Kondo effect in asymmetric Josephson couplings through a quantum dot 
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Asymmetry in the Josephson couplings between two superconductors through a quantum dot 
is studied based on a single impurity Anderson model using the numerical renormalization group 
(NRG). Specifically, we examine how the difference between the couplings Tl and Tr affects the 
ground state, which is known to show a quantum phase transition between a nonmagnetic singlet 
and a magnetic doublet depending on the various parameters; the Coulomb interaction U , onsite 
potential e^, level width caused by the hybridization, and superconducting gap Aj, = lAj^le*^" 
for the leads on the left and right [v = L, R). Our results show that whether the local moment is 
fully screened or not depends substantially on the asymmetry in the couplings Tl F_b. It tends to 
make the singlet ground state stable, while the amplitude |A^| and phase difference (j> = Or — 9l of 
the superconducting gaps tend to suppress the screening. We also discuss some general symmetry 
properties of the system and their relation to the current conservation. 

PACS numbers: 72.10.-d, 72.10.Bg, 73.40.-c 



I. INTRODUCTION 

The Kondo effect in superconductors has been stud- 
ied over three decades ^ii^i^i^ It was shown in early years 
that the competition between the superconducting (SC) 
gap and Kondo energy scale Tk determines the low- 
temperature properties. The SC energy gap Asc dis- 
turbs the conduction-electron screening of the local mo- 
ment. The ground state becomes a magnetic doublet for 
^sc ^ Tk, while the ground state is still a nonmag- 
netic singlet for Asc ^ ?k- These aspects of the Kondo 
physics in superconductors have been re-examined pre- 
cisely by efficient numerical methods such as the quan- 
tum Monte Carlo (QMC)^ and numerical renormaliza- 
tion group (NRG) approaches i^'^'^ 

In quantum dots embedded in a Josephson junc- 
tion, new interesting features have been predicted to 
be observediSiiiii^ and some experiments have been 
reportedji^ii^ii^ii^ Particularly, the screening of the lo- 
cal moment in the quantum dot is affected not only by 
the size of the two SC gaps but also by the Josephson 
phase (f) between the two SC leads on the left and right. 
The phase difference 4> also induces the Josephson cur- 
rent flowing through the dot. It can drive a quantum 
phase transition (QPT) between the singlet and doublet 
ground states, and at the critical point the direction of 
the current changes discontinuously. 

These properties of the SC-dot-SC systems have been 
studied with various theoretical approaches; such as the 
noncrossing approximation, slave-boson mean-field 
theory ; -'^^1^° perturbation theory in the Coulomb inter- 
action QMC^i23 and NRC^i^s So far, however, 
most of the calculations have been carried out assum- 
ing a highly symmetric condition as |Al| — \^r\ and 
Tl = rH^a8a9^>2i^^ jjere^ ^ |A,|e'«- is the 

SC gap of the leads on the left and right {v — L, R), 
and is the bare level width due to the hybridization. 
Therefore, the feature of the QPT far from the symmet- 



ric condition is one of the theoretical issues to be clarified 
with some accurate approaches. 

In a previous work, we have studied the ground-state 
properties in the case of |Al| |Afl|. In real quantum 
dots, it corresponds to the situation that the two su- 
perconductors connected to the dot are not identicali^ 
Specifically, we have considered the limit where one of the 
two SC gaps, Al, is much larger than the other one A^j. 
In the limit of |Al| — > oo, the left lead has been shown 
to be separated from the rest of the system, which con- 
sists of the dot and right lead. It just leaves a static SC 
pair potential, which is given by A^ = e*^^ , in the 
dot. The pair potential appears because a SC proximity 
from the left lead is allowed even in the large gap limit, 
while the quasi-particle tunneling into the left lead is pro- 
hibited. This situation is described by a single-channel 
Anderson model with an extra SC pair potential A^ in 
the impurity site. The simplified model possesses essen- 
tial ingredients of the Kondo and Josephson effects, and 
the phase difference (j) = 9fj — 9l can be defined between 
Ad and A^,;. We have studied the features of QPT ap- 
plying the NRG to this model, and have confirmed, for 
example, that the Andreev bound state emerging close to 
the Fermi level causes a re-entrant QPT near half-filling 
for (/) ~ TT and Tl ^ Tr^ 

Another aspect that seems to be important is how the 
asymmetry in the couplings and F^ affects QPT. 
That is the main subject of the present work. The SC 
proximity effect depends strongly on the environment 
surrounding the dot, and the competition between the 
Kondo and proximity effects determine the low-energy 
properties. In this paper we present the results of the 
QPT phase diagram obtained for the asymmetric cou- 
plings Tl ^ Ffl. The NRG calculations have been car- 
ried out mainly using the two-channel model, assuming 
the two gaps are equal in size |Al| = |Afl|. Our re- 
sults show that the asymmetry in the couplings tends to 
enhance the screening of the local moment. 



The model and some symmetry properties of the An- 
derson impurity in the Josephson junctions are described 
in SecHJ based on a rotational symmetry in the pseudo- 
spin space. The NRG approach and the results of the 
phase diagram for QPT obtained in some typical param- 
eter regions are discussed in Sec. IIIII We also discuss in 
Sec. llVl the roles of the asymmetry in two solvable cases: 
noninteracting case U = 0, and a superconducting ver- 
sion of the atomic limit for |Al| — > oo and \A.ji\ — s- oo 
which can also be considered as one of the fixed points 
of NRG. A summary is given in Sec. |Vl We also de- 
scribe some general properties deduced from the current 
conservation in the appendix. 



II. MODEL FOR A SC-DOT-SC JUNCTION 

We describe in this section the model and some of the 
symmetry properties. Specifically, at half-filling for the 
Josephson junctions with = or tt, the system has a 
global U(l) symmetry in a pseudo-spin space. It enables 
us to map the SC leads onto the normal conductors with 
a staggered potential, as shown by Satori et al^ 



A. General formulation 

We start with a single Anderson impurity connected to 
two SC leads on the left (L) and right (i?), as illustrated 
in Fig. [TJ The Hamiltonian is given by 
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Here, S,d = + U/2, with ed the impurity level, and U 
is the onsite Coulomb interaction. The operator dj. cre- 
ates an electron with spin a at the impurity site, and 
Ud = X^CT^CT^cr- The other operator c^q^ annihilates a 
conduction electron that couples directly to the dot via 
the tunneling matrix element v^^. The electrons in the 
leads can be described by the Hamiltonian of a tight- 
binding form without loss of generalities.^^ The hop- 
ping matrix element „ determines a continuous one- 
particle spectrum ekv and eigenfunction (fj^^in), by which 
the hybridization strength can be written in the form 



A/. 



Vl^ Vr 



FIG. 1: Schematic picture of an Anderson impurity con- 
nected to two SC leads with the gaps Al and A_b, where vl 
and Vr are the tunneUng matrix elements. 



The symmetry properties of the system can be clearly 
seen, in the Nambu representation using the spinors 



(-l)-'^c 



Then, eqs. (H))-® can be rewritten in the form 
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Here, Aj, and id are defined by 
A* 



id 



(6) 
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and Tk for fc = 1, 2, 3 is the Pauli matrix. The Coulomb 
interaction Ti.^ is isotropic in the pseudo-spin space, and 
the expression given in eq. ([7]) , which is identical to that 
in eq. ^ , explicitly shows the rotational symmetry. Fur- 
thermore, the current through the dot can be expressed 
in the form, 

^« = - * ^ ( ^R.o ^3 V-d - ^3 rpjifi ) - (12) 

JL = -^'^{rplT,tP^,-xPl„T3rP^). (13) 

The current is conserved such that dn^/dt + Jj^ — J]^ — 0, 
as described also in the appendix [Al 

The rotational symmetry in the pseudo-spin space are 
described by the operator 



id 



oo 

9 E E^^^«^^ 



(14) 



i'=L,Rn=0 



The z component corresponds to the half of the total 
charge Q = 21^, and Ix [ly) represents the real (imagi- 
nary) part of the SC pair potentials. These three com- 
ponents satisfy the same commutation relations as that 



3 



IA/.I 



(p=0 



Q 



(a) 

IArI 



'c/or 



(p=7: 

n 



(b) 

IArI 



□ 



'c/or 



0. The Horizontal line 



FIG. 2: Potential profiles for the quasi-particles ^i^,na, given 
in eqs. (|16|l - (|20p . for the Josephson junctions with (a) 4> = 0, 
and (b) tf) — n, ai half-filling 
corresponds to the Fermi level. 



of the angular momentums, 



[Iy,h]^tl,, [h,l,]=ily. (15) 



In the pseudo-spin space, the onsite potential can be 
interpreted as a local external field along the z axis, and 
the SC gaps can be regarded as a staggered field in 
the x-y plane. Therefore, if = and Ai, = 0, which 
corresponds to normal leads in the particle-hole symmet- 
ric case, the system has a full rotational symmetry in the 
pseudo-spin space as well as the rotational symmetry in 
the real spin space, because [7^, 7i] = for k = x, y, z. 

Another typical case is that of the Josephson junction 
Ay 7^ with the particle-hole symmetry = 0. In 
this case the z component of the external field vanishes. 
Therefore, if we take the quantization axis to be in the 
y-direction by carrying out a spinor rotation of 7r/2 with 
respect to the x axis, the Hamiltonian Ti. takes a real 
symmetric form. Then, with the new axis, the current 
operators are described by a Pauli matrix T2 instead of 
T3 in eqs. (dm and p^ . 



B. Particle-Hole symmetric case with real gaps 

The two SC gaps A^ and A/j are essentially real for 
the Josephson junctions with = or vr. In this case, 
if the system also has the particle-hole symmetry = 
0, the model can be mapped onto a normal asymmetric 
Anderson model with a staggered potentiali^ It is owing 
to a conservation associated with [Ix,H] = 0, caused 
by a U(l) symmetry along the x axis. It enables us to 
transform into the total number of the quasi-particles 
as shown in eq. (j2ip . This feature can be seen explicitly 



if we take the quantization axis to be in the x direction^ 
carrying out a spinor rotation of 7r/2 with respect to the 
y axis by a unitary transform e'~'^^ as 

1 1 

(16) 

Then, the Hamiltonian can be expressed in the form 
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(17) 
(18) 



v=L.Rn=Q a 

(19) 

= E E(-l)"+^Ai^7i,„.7..„.-l). (20) 

u=L,Rn=a \ o- / 

Note that H^^ = 0, because of the particle-hole symme- 
try; = 0. Therefore, the energy spectrum is deter- 
mined by the quasi-particles, the total number of which 
is conserved 

24 = + E E ( E^i,-^-,- - H ' (21) 



(22) 



Particularly, the local charge q_i corresponds to the SC 
pair correlation penetrated into the impurity site. The 
profile of the potential that the quasi-particles feel is 
illustrated in Fig. [2l For the Josephson junction with 
(f) — the potential still shows a regular staggered period 
around the dot, while for </> = tt the steps of the poten- 
tial changes like that in a defect, or dislocation. The 
potential profile suggests that the SC proximity effect, 
which causes a finite anomalous average {q_i), is larger 
for = than that in the case of tt. The penetration 
of the SC correlations into the impurity site depends on 
the tunneling strength F^ and the size of the SC gaps 
A^. Specifically, for (j) — tt and F^ = F^, the impurity 
level evolves to the Andreev bound state emerging just 
on the Fermi level situated in the middle of the energy 
gap. In this case, if also the size of the two gaps are equal 



A^|, the SC correlations cancel out (q^i) = in 



the impurity site. It also holds in the particle-hole asym- 
metric case ^ as shown in the appendix [X] based on 
the current conservation. 



C. Particle-Hole asymmetric case for 



There is a different symmetry for a 7r-junction (f> = n, in 



the case of F 



R 



F^ and 
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It holds without 
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the particle-hole symmetry, and reveals, for instance, in 
the Nambu form of the Green's function G given in eq. 
(jA12p in the appendix|Al The underlying symmetry links 
to the conservation of a channel charge Qch defined by 



bLb. 



n—0 a 



V2 



V2 



(23) 



The commutation relation, [Qch, = 0, can be con- 
firmed explicitly if we express eqs. dS])-® in terms of 
bonding and antibonding orbitals for vl = vr, and tak- 
ing t^^n to be channel independent. The superconducting 
part, eq. ([5]), conserves the difference between the total 
number of the particles with the even symmetry and that 
of the odd ones. 
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|Al| = |Ar| (=A), 8d = -U/2 
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FIG. 3: Phase diagram of the ground state for = at half- 
filhng = —(7/2 is plotted for several values of A: (□) 1.5, 
(•) 2.0, (A) 2.5, (o) 3.0, and (A) 6.0, where Vl = Vr (= T), 
and |Ai| = |Ah| (= A). 



Here, A is real and defined such that A^ — — — 
A. The Hamiltonian can be simplified carrying out the 
particle- hole transformation for 



Then, the two leads consist a normal ladder 



1+1 ?,t 



n—0 (J 

oo 



(25) 

(26) 
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n—0 (7 



and one of the corners for is connected to the dot 

as Ht — X^CT V2^('*0(t'^<t + ^i'^ocr)- The dot part of the 
Hamiltonian, TiJ] + 7i^, keeps the original form given 
in eq. (O. Thus, it is obvious that the SC correlation 
vanishes, (dlrfj) = 0, at the impurity site in the present 
situation even for the particle-hole asymmetric case ^ 
0. 

The ladder described by eqs. and (P7)) gives the 
density of states that shows a square root divergence at 
the edges of the gap e = ±A. Therefore, the properties 
of the TT-junction in the symmetric gaps and couplings 
must be essentially the same as those of the magnetic 
impurity in a band insulator, which has been studied 
using mainly a constant density of states with the gap 
around the Fermi levelj^i^liS^ 



III. GROUND-STATE PROPERTIES 

A. NRG approach 

We study the ground-state properties using the NRG. 
In order to examine the effects caused by the asymmetry 



in the tunneling matrix elements, we consider the Joseph- 
son junction with |Al| = \^r\ (= A), where the two gaps 
are equal in size. To be specific, we assume that r^(e) is 
a constant inside the conduction band —D < e < D, and 
use the standard logarithmic discretizationi^ Then, the 
hopping and tunneling matrix elements are expressed in 
the form 

I + IM (l-l/A"«)A-^' 

2 ^1 - l/A2"+l^l - l/A2«+3 

I 1 + I/A 



2Dr^AA 



2 1 - 1/A 



log A. (29) 



The factor A\ is required for reproducing the original 
model correctly in the continuum limit A — > li^i^ 
We have carried out the successive diagonalizations of 
NRG keeping typically the lowest 500 energy states for 
constructing the Hilbert space with extra orbitals from 
next step. The low-energy properties are determined 
by the ratios T„/A, U/A, e^/A and the phase differ- 
ence (j>. In the calculations we have taken the gap to 
be A/D = 1.0 X 10^^: this ratio itself does not affect 
the low-energy properties, provided A ^ D. We have 
used the eigenvalue of the pseudo spin as a quantum 
number to classify the Hilbert space in the particle-hole 
symmetric case for (f) = and tt. These two limits of (j) 
provide the upper and lower bounds in the ground-state 
phase diagrams for intermediate (j). 

The model for the Josephson junction contains two 
conduction bands, and the accuracy of the NRG results 
becomes somewhat worse than that for a single-channel 
model. For this reason, we have mainly used a rather 
large value A = 6.0 for the discretization parameter to 
get an early convergence in the NRG steps. To check the 
validity of it, we have examined the phase diagram of the 
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FIG. 4: Number of electrons in the dot (rid) as a function of 
ea for U = 5.0A, |rt| = |rfl| {= T), and |Ai| = \Ar\ {= A). 
The ground state is a singlet for (•) F — 0.35A, while for (A) 
F = 0.65A it changes to a doublet for -4.2 < ed < -0.8. The 
discretization parameter is taken to be A = 2.0. 



FIG. 5: Phase diagram of the ground state at half- filling ed = 
— 17/2 for asymmetric couplings Fl ~ 1.44F_r with (•)</!> = 
and (a) = tt. The dashed lines show the phase boundaries 
for the symmetric couplings Fl = F_r with (o) (p = and (A) 
(j> = Ti. The SC gaps are taken to be |A_l| = 1A_r| (= A). 



ground state changing the value of A (= 1.5, 2.0, 2.5, 3.0 
and 6.0): we have preserved the lowest 2000 states for 
A = 1.5. The calculation has been carried out for 
(/) = where the model can be reduced to a single-channel 
model, and the results obtained at half- filling = — C//2 
for the symmetric couplings T l = ^ r (= T) are shown 
in Fig. [H We see that the results obtained at A = 6.0 
already agree well with those obtained for small A's. It 
indicates that the level crossing of the ground and first 
excited states is not so sensitive to A, and it supports the 
validity of our results presented in the following. The fig- 
ure [3] itself shows that the ground state is a nonmagnetic 
singlet for large F or small [/, and at the critical point it 
changes discontinuously to a magnetic doublet for small 
F or large 

In order to show an example of the behavior of ob- 
servables at the critical point, we have also calculated 
the average number of electrons (n^) in the dot. The re- 
sult is shown in Fig. [4] as a function of for U — 5. OA, 
= 0, and Tl = Vr. In this parameter set, the ground 
state for F = 0.65A is a singlet for all values of e^, and 
thus {rid) decreases monotonically with increasing e^- In 
contrast, for F = 0.35A, the ground state is a doublet 
for —4.2 < Ed ^ —0.8, and (n^) jumps at the critical 
points. Other observable quantities, such as the Joseph- 
son current, will also show similar discontinuous changes, 
because the QPT between the singlet and double ground 
states is the first order transition occurring as a result of 
a level crossing of the Andreev bound state emerging in 
the SC gap in the energy spectrum. 



B. Asymmetric couplings Vl ^ F^ at half- filling 

How does the asymmetry in the couplings F^ and F^^ 
affect the ground-state properties? We first of all con- 
sider the particle-hole symmetric case ~ — C//2. In this 
case, the dot is occupied by a single electron (n^) = 1 on 
average, and magnetic correlations are enhanced maxi- 
mally. In Fig. [51 the phase diagram of the ground state 
for asymmetric couplings F^ = 1.44F/J is plotted for the 
Josephson junctions with (•)(/) = and (a) <j) = tt. For 
comparison, the results obtained for the symmetric cou- 
plings Tl = Ffl are shown for (o) = and (A) (f> = tt 
with the dashed lines. The phase boundary for interme- 
diate values of 0, namely for < 4> < n, will appear in 
between the line for = and that for (j> = ir. In the 
figure, we see a general trend that the phase difference (j) 
between the SC pair potentials enlarges the region for a 
doublet ground state. It means that the Josephson ef- 
fects due to a finite 4> disturbs the screening of the local 
moment, and the QPT will be driven by if the other 
parameters U/ A and Fj^/A bring the system into the pa- 
rameter region enclosed by the two lines. Note that for 
the data obtained for Fl = IAATr, the sum of the hy- 
bridization energies F^ + Tr is larger than that for the 
symmetric couplings. It also favors the singlet ground 
state. 

The results shown in Fig. [5] suggest, however, that the 
superconducting proximity effects on the impurity site 
depends substantially on the asymmetry in the Josephson 
couplings, particularly for the 7r-junction. The vertical 
phase boundary at = tt for F^ = F^ is caused by an 
Andreev bound state. It emerges just on the Fermi level 
at half-filling, as we see also in the potential profile shown 
in Fig. [21(b). Therefore, along the vertical line the ground 
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FIG. 6: Phase diagram of the ground state away from half- 
filhng at ed = — 2.5A is plotted for asymmetric couplings Fl ~ 
IAAFh for {•) <j> = and (a) (f) — n. Corresponding results 
for the symmetric couplings Fl = Tr are also plotted for 
(o) (j) — and (A) (p = n. Particle-hole symmetry holds at 
U = 5.0A. 
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FIG. 7: Phase diagram of the ground state away from half- 
filling at [/ = 5. OA for asymmetric couplings F_l = 1.44F_r for 
(•) (f} = 0, (a) (f) = TV, and the symmetric couplings Fl ~ F_r 
for (o) (p — 0, (A) (p — n. Particle-hole symmetry holds at 
ed = -2.5A. 



state has 4-fold degeneracy. An infinitesimal U lifts the 
degeneracy when the two SC gaps are identical |Ai| = 
and the ground state becomes a doublet in the 
whole region of C/ > 0. We will discuss this point more 
in detail in Sec. IIVI The asymmetry in the couplings 
Tl j^Tb. breaks this situation. Then, the Andreev bound 
state emerges away from the Fermi level, and it makes 
the ground state to be a singlet for small U. 



C. Away from half- filling 

We also consider the ground state away from half-filling 
7^ U/2. The phase diagram obtained for a fixed value 
of the onsite energy — — 2.5A is shown in Fig.[6l where 
the particle-hole symmetry holds along the dashed line 
at J7 = 5. OA. The calculations have been carried out 
for asymmetric couplings Fl = 1.44Ffl for {•) (j) — 
and (a) (j) — ""i E^nd also for the symmetric couplings 
Tl = for (o) = and (A) (j) — n. As mentioned, 
the boundary for < < tt will appear between the line 
for (j) — and that for cf) — tt. In Fig.[6l all the bound- 
aries start from the point U — 2.5A in the horizontal 
axis. This is because in the atomic limit = the 
level crossing between the singlet and doublet states oc- 
curs at 2ed + U — ed, when the energy of single-particle 
state and that of two-particle state coincide. We also 
see in this figure that the phase boundary for (j) — tt 
is changed substantially by the asymmetry in the cou- 
plings, and the parameter region corresponding to the 
singlet ground state spreads in an upper-right region for 
Tij ^ 1.4A and U > 5. OA. For the symmetric couplings, 
the phase boundary at = tt approaches asymptotically 



to the vertical line at U — 5. OA. This behavior is caused 
by the Andreev state approaching to the Fermi energy. 

In quantum dots is a tunable parameter that can 
be controlled by the gate voltage. The dependence 
of the phase diagram of the ground state for F^ = Tr 
has been obtained first with the slave-boson mean-field 
approximation by Rozhkov and ArovasJ^ We have ex- 
amined the effects of the asymmetry. The NRG results 
are shown in Fig. [7] for (solid line) the asymmetric cou- 
plings Fi — 1.44Fi^, and for (dashed line) symmetric 
couplings Tl — Tfl. Here, the Coulomb interaction is 
fixed to he U = 5. OA. The ground state is a singlet for 
ed + U < and > 0, because the impurity site is al- 
most fully filled or almost empty in these two regions. We 
see again in Fig. [7] that the asymmetry in the couplings 
Tl 7^ enlarges the parameter space of the singlet 
ground state. The phase boundary for </) = tt in the sym- 
metric couplings Fi = Ffl shows a diverging behavior at 
ed — — 2.5A, which is owing to the Andreev bound state 
being situated just on the Fermi level at half-filling. 

IV. DISCUSSIONS 

The phase transition between the singlet and doublet 
ground states occurs at a level crossing of the lowest and 
the first many-body excited states, both of which are dis- 
crete and are separated from the continuous excitation 
spectrum at e > min(|Aj^|, |Al|) above the energy scale 
of the SC gap. Therefore, a few discrete states emerging 
inside the gap region determine the properties at low en- 
ergy scales. In this section, to obtain some insights into 
how the asymmetry around the dot affects the ground 
states, we consider two simple limits where the model 
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FIG. 8: Ground state for (f> — tt, Fr = Fl, and S,d — 
in the presence of an infinitesimal Coulomb interaction U, 
which lifts the degeneracy caused by the zero-energy Andreev 
state. If ere, x = [Afl|/|Az,|, represents the asymmetry in 
the SC gaps, and \Al\ is changed in the region |Ai| > \Ar\ 
keeping |Ajf,| unchanged. The critical value of F at a; = is 
F„ ~ I.77|Afl|. 




X 

FIG. 9: Bound states for U — exist in the region upside 
of each of the lines for (p — 0, n/6, 7r/4, and n/3. For n/2 < 
(j) < n, the bound states always exist. Here, x — |A_r|/|Al|, 
and |A_t| is changed in the region |Al| > |A_r| keeping \Ar\ 
unchanged. The lines are determined by eq. p2p . which does 
not depend on F_b nor 



can be solved. 

A. Noninteracting limit Ti.^ = 

We consider here how the asymmetry in the two SC 
gaps affects the Andreev state in the noninteracting case 
— 0. The density of states of each SC lead has the 
square root divergence at the edge e = ±|Aj^|, and the 
formation of the bomid state depends strongly on these 
divergences. The edges are separated from each other for 
|A^| 7^ |A^|, and they merge for the symmetric gaps to 
give a linear divergence for < < tt. 

The energy of the Andreev state inside the gap is deter- 
mined by an equation, det {G''(e)} = 0. The Green's 
fmiction is defined in the appendix [XI The explicit form 

of the determinant F{e) = det{G'^{e)} is given, re- 
placing the imaginary frequency iun in eq. (|A3[) by the 
real one e for |e| < mindAflj, |Al|), as 

VWIf^ y/\An\^-e^ 
2r^r;,.(6^-|Az.||A;,.| cos^) 
v/|Ai|2-eVIAflP-e2 

It takes a negative value, or vanishes, at e = 0, 

(31) 

Therefore, the bound state emerges if F{e) increases 
with |e| to reach zero. For 7r/2 < < tt, the solu- 



tion always exists because in this range of the phase 
-^(e) is an increasing function and diverges positively at 
|e| = min(|Aij|, |Al|) with the square root dependence. 
Specifically, the bound state appears just on the Fermi 
level e = 0, when the right hand side of eq. (pij) equals 
zero. It happens for = 0, = tt, and T^ = F^, pro- 
vided |A^| and |A^| are nonzero. In this case the sin- 
glet and doublet ground states are degenerate because 
the energy does not depend on the way the bound state 
is occupied. The asymmetry in the couplings F^ and 
F^ affects crucially the ground-state property for the tt- 
junction. The Coulomb interaction lifts the degeneracy, 
and it can be studied with a simple perturbation the- 
ory with respect to an infinitesimal U ^ The results are 
shown in Fig. [51 assuming that IAt^J = min(|A7?|, |Ai|). 
The way how the degeneracy is lifted depends quantita- 
tively on the asymmetry in the SC gaps x = |A^|/|Aj;^|. 
Specifically, the curve for the phase boundary reaches 
at a; = 1, which corresponds to the symmetric gaps 
|A^| = |A^|. It means that the ground state is a dou- 
blet for any values of F at a; = 1. It explains the behavior 
of the vertical phase boundary for the 7r-junction in Fig. 
[5l and in this case the SC correlation in the impurity 
site vanishes (rfjdj) = also for finite U as shown in the 
appendix [^ 

Another example seen in the case of asymmetric gaps 
|A^| 7^ |A^| is that the bound state disappears in a 
finite parameter region. It is caused by the difference in 
the asymptotic behavior F(e) at |e| |A^| — 0+, where 
we assume that |Aj^| < |A^|. For instance, in the limit 
of |A^| oo, the determinant F{e) does not have a 
zero point for |A^| < F^ cos(/>!^ This is because the SC 
correlation of the order F^, which penetrates from left 
lead into the dot, pushed the bound state away towards 
the continuum outside of the gap region of the width 



|Aj^|. The bound state can remain inside the gap for 
|A^| > Tj^coscj). The same situation also happens for 
finite |A^|. The existence of the bound state can be 
deduced from the asymptotic behavior of the determinant 
at the edge of the gap region. The bound states appears 
if F{e) +00 for |e| |A^| — 0+, and this condition 
can be expressed in the form, for < a; < 1 with x — 
|Ai^|/|Ai|, as 

|A^| cos(f,~x 



The solution does not exist when the both sides of eq. 
([5^ are equal, since F{e) takes a negative finite value for 
|e| \^r\ — 0"*" at the critical condition. The param- 
eter region where the bound sates do not exist shrinks 
as 4> increases as shown in Fig. [S] In contrast, for the 
symmetric gaps |A^| = |A^|, the determinant shows a 
positive linear divergence for < < tt at the edge of 
the gap, and the bound state always remains inside the 
gap region. 

B. Large gap limit |An|, |Al| 00 
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FIG. 10: Equi-energy surface of Andreev state Ea in a pa- 
rameter space consisting of Tn, Fl, and in the limit of 
|Ai{|, |A_l| — > 00. The phase </> changes the eccentricity of the 
ellipse in the horizontal plane, and the picture shows an ex- 
ample for a small (j> (7^ 0). The central region between the two 
white longitudinal lines, which represent F_r — and Tl ~ 0, 
is the physical parameter space for F^ > 0. The ground state 
is a doublet inside the surface which is determined such that 
Ea = U /2, and is a singlet for the outside. 



When both of the two gaps |Ajj| and |Aj^| are larger 
than the other parameters, the model can be mapped 
onto a simple single site problem by taking a limit of 
|A^| 00 and |Aj^| ^ ooi^ It can be regarded as a 
superconducting analogue of an atomic limit, and can 
be solved analytically as summarized in the appendix 
IB] The properties in the symmetric couplings F^ — F^ 
have been studied by Vecino et al^^ In the following, we 
discuss mainly the role of the asymmetry in the couplings 

In the case of a normal transport through a quantum 
dot the equation of motion method,^ which is essentially 
the atomic limit in the Hubbard I approximation^^i has 
been applied to the properties at high energy scale. It 
has a limitation, however, at low temperatures because 
it does not describe properly the Fermi liquid proper- 
ties due to the gapless low-energy excitations. Never- 
theless, in the case of the dot connected to supercon- 
ducting reservoirs, the low-energy excitations have the 
gap, and it gives a different meaning to the atomic limit. 
The single site degrees of freedom can be identified as 
the Andreev bound state remaining inside the gap. It 
does describe the lowest few many-body excited states 
inside the gap properly, and can also be considered as 
one the fixed points of NRG. Although the virtual pro- 
cesses to the continuous excited states above the SC en- 
ergy gap may also be important for the correlation func- 
tions, the main features of the contributions of the An- 
dreev state to the low-energy properties, especially the 
bound-state energy itself, can be captured in the atomic 
limit of \Al\,\Ar\ 00. 

The effective Hamiltonian for the single site model is 
given in eq. (|B4p in the appendix [B) In the large gap 
limit, the role of the SC leads can be replaced by the SC 



pair potential at the impurity site A^ — F^e*^^-|-F^e'^^, 
which is caused by the proximity effects. The energy of 
the Andreev state Ea = + l^dp given in (|B5p can 
be rewritten in the form, 

EA^^ea + i^R + r J' cos2 ^ + (F^ T^f sin^ t . 

(33) 

Therefore, the equi-energy surface of Ea can be described 
as an ellipsoid in a three-dimensional space consisting of 

and F^ ±F^ as shown in Fig.[TOl The Josephson phase 
(j) changes the aspect ratio, or eccentricity, of the ellipsoid. 
Specifically, for (j) ir^ the gradient along the F^ — F^ 
axis becomes large, and it implies that the asymmetry 
in the coupling plays an important role on the properties 
near the 7r-junction. 

The ground state is a doublet for Ea < J7/2, and it is a 
singlet for Ea > U /2. Therefore, the Coulomb interacion 
determines the size of an ellipsoidal such that 

Ea = U/2 . (34) 

The inside corresponds to a doublet ground state, and the 
outside is a singlet. The equations ([55]) and explain 
the features of the phase diagrams, shown in Figs. EHTl 
for small parameter values, at max(F^, F^, [/, \S,d\) ^ A. 
Particularly, Fig. [7| can be related to a cross section of 
the elliposoid for a fixed U/2. For instance, for (p = tt, 
the phase diagram in the atomic limit is determined by 
frf + (F^ - F^)^ = and thus the phase boundary 

becomes two open sheets at = ±11/ 2 for the symmetric 
couplings Fjj. = F^, while the asymmetry F^ ^ F^ closes 
the sheets to form an elliptic cylinder. As another simple 
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example, we consider the symmetric couplings = 
(= F), for which the boundary is given by 4F^ cos^ ^ = 
[7^/4 — This equation corresponds to the dashed lines 
shown in Figs. [H] and [5] for small F. 

In the atomic limit the first and second many-body ex- 
citation energies Ei and E2, measured from the ground- 
state energy, are given by 



singlet GS: 
doublet GS: 



Ei=Ea- U/2 
E2 = 2Ea 

E, = U/2 - Ea 
E2 = U/2 + Ea 



for Ea > 



for Ea < 



U 
2" 

u 



(35) 
(36) 



Two discrete many-body excited states can appear in the 
energy gap even for finite SC gaps, for which the ener- 
gies can be calculated with NRG. Then, eqs. ([35|) and 
P6p could be used to determine the renormalized val- 
ues of U and Ea, which asymptotically coincide with the 
bare ones in the limit of \Aj^\, |A^| — s- 00. The renor- 
malization is caused by the contribution of the cotinuous 
states above the energy gap. The analysis of the fixed 
point along these lines will be discussed elsewhere. 



V. SUMMARY 

We have studied effects of the asymmetry in the 
Josephson couplings F^ 7^ Tr on the ground state of 
a quantum dot embedded between two superconductors. 
Specifically, we have obtained the phase diagram of the 
ground state in some parameter regions. The results 
show that the asymmetry in the couplings enhances the 
screening of the local moment in the dot, particularly for 
the TT junction. For ~ tt, the Andreev bound state 
emerges near the Fermi level in the symmetric couplings 
Tl = Ffl;, and it favors the magnetic doublet ground state 
for U > 0. The asymmetry in the couplings moves the 
Andreev bound state away from the Fermi level, and en- 
hances the screening of the local moment. The boundary 
between the singlet and doublet ground states for inter- 
mediate phase < < TT must be in between the curves 
for (/) = and that for (p = n. 

We have also described some symmetry properties 
of the Anderson impurity in the Josephson junction in 
terms of the spinor rotation in the pseudo-spin and the 
current conservation. Specifically, the similarity between 
the TT-junction and a magnetic impurity in an insulating 
host has been dicussed based on a symmetry property. 
Furthermore, the role of the asymmetry in the couplings 
and gaps has been discussed in the noninteracting case, 
and also in the hmit of |A^|,|A^| 00. Our results 
suggest that in real quantum dots the ground-state prop- 
erties depends in various ways on the local conditions 
surrounding the dot. 
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APPENDIX A: CURRENT CONSERVATION 

The imaginary-time Green's function for the Anderson 
impurity in the SC hosts can be described in the Nambu 
representation as 



G(iLu„ 



P 



'{TrdJT)d^^ {TrdjT)d^) 
{Trdl{T)d\) (Trd\{T)d^) 



(Al) 



where tj„ = (2?i + l)7r//3 for n = 0,±1,±2,. .., /3 = 1/T, 
and T the temperature. The Dyson equation is given by 

{G{lLJn)r' = {G°(iC^„)}"' - ^{iuJn) , (A2) 

where S(ia'„) is the self-energy caused by TC^ , and 

{G (iLUn)} = lUJni ~ idTs + F^ 



(A3) 

The currents through the dot is defined based on the 
Heisenberg equation, dud/dt = {^i/fi) [nd, Ti] , as 



ot 



Jr 



Jr. = 



(A4) 
(A5) 

(A6) 



Here, Jl is the current fiowing into the dot from the left, 
Jr is the current from the dot to the right. The thermal 
average of these currents can be expressed in terms of the 
Green's function 



(Jl) 



2e 



E 



-iT, 



Tr 



Tr 



T3Aj^G(«w„)J, (A7) 
t^AlG{iuj^)]. (A8) 



Here, the trace is taken over for the 2x2 matrices for the 
pseudo spin space. The current conservation ( J^,) = ( Jj?) 
indicates that there is one relation between the off diago- 
nal elements of G(ia;„). For instance, the self-energy ma- 
trix can be described by three functions ao(iw„), a3(ia'„). 
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and a-^g.y^piiujn^iujyn) ■ 



\ Sj.(iw„) I 



(A9) 



(AlO) 



pA 



M4 M2,P3 



The kernel a'^^.-^p{iuimi^m) for the ofF-diagonal part of 
the self-energy Sx(ia;„) can be shown, using a Ward- 
Takahashi identity;^ to be expressed in the form, 

x{GilLUm)T3} {CiliUra)} ■ (All) 

(iuj4, iuji ; ia;2, ^^3) is the vertex function 
in the Nambu representation, and [j = 1,2,3,4) rep- 
resents the component of the pseudo spins. From these 
properties, it is deduced that, if the both couplings and 
gaps are symmetric — F^ (= F) and \/^r\ — |Al| 
(= A) for the 7r-junction = tt, the off diagonal part 
of the the self-energy vanishes at the impurity site. It 
simplifies the Green's function as 



Here, F 



P4P1 ;P2P3 



{G{iuj,-,)}-^ 



1 



2F 




^^l + A2 
1 + a3(icj„)^ T3. 

Therefore, in this particular case the SC correlation can- 
cels out in the impurity site = 0, as can be deduced 
from eq. (|Aip . Furthermore, if the system also has the 
particle-hole symmetry = additionally, the Green's 
function is simply proportional to the unit matrix 1. 

APPENDIX B: LIMIT OF \Ar\, Al ^ 00 

In the limit of |Al| — > 00 and |A;j| — > 00, the model 
can be mapped onto a simple single-site problem fSii^i^^ 
which can be solved analytically. In this limit, the Dyson 
equation eq. (|A2[) can be rewritten in the form 



id Ad 



[iuJn) 



A* 



A, 



F^e'^-. 



(Bl) 
(B2) 



The SC correlation penetrates into the dot from the SC 
leads via F^ and F^^. Therefore, the amplitude depends 
on the Josephson phase (j) = Or — 9l, 



\Ad\' = {Tj 



rr)'-4F^FiSin2| 



(B3) 
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FIG. 11; Statistical weight p_ — p+ as a function of Ea for 
r/17 = 0.0, 0.1, 0.2, and 0.5. 



The corresponding single-site Hamiltonian is given by 

2C/ ^ ^ 



(B4) 



The noninteracting part determines the energy scale of 
the Andreev bound state Ea = \J(\ + lA^p. The 
ground state is a singlet for > U/2, while it is a 
doublet for E^ < U/2. The boundary between these two 
ground states is determined by the equation U/2 = Ea- 
Since the bound-state energy Ea depends on (j) — 6ji — 6l 
as 



Ea = \jed + (r^ + ^l) - "^^B^L sin^ f , (B5) 

the phase transition occurs in between < < tt, if the 
Coulomb interaction is in the region Eatt < U/2 < Eao- 
Here, Eao^tt is the value of Ea at = 0, tt. 

The Green's function can be obtained analytically 



> ■c^a + (i + .)ife)HS 



P± = 



2 + e~'^(t+SA) +e-''(¥--^'*) ' 
The equal time correlations, G(r)|^^_Q+, is given by 



(B6) 
(B7) 



{dt dj) (d^ d^) 



= -1 
2 



2Ea 



H 



d ■ 



(B8) 



Note that in the noninteracting case U — 0, the statistical 
weight p_ — is replaced by 1 — 2/(i?^) with /(e) = 

[e^' + iy ■ InFig.ini the weight p_ — is plotted as a 
function of Ea/U for several temperatures. Specifically 
at zero temperature, p_ — I and = 0, for Ea > U/2 
in the singlet ground state. It takes the value — = 
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1/2 for Ea < U/2 in the doublet ground state. Therefore, 
at T = 0, the self-energy can be expressed in the form 



U 
2Ea 



U 



2 _ r2 



(B9) 



The Josephson current can be calculated substituting 
eq. (HU into eq. (|T7)) or 



(^)-^^4r,r,sinc 



(BIO) 



where Eaq = + (^-R + ^h)^ is the maximum value 
of the bound-state energy, and T = ^ j^j E\^ cor- 
responds to a normal transmission probability for non- 
interacting electrons. The contribution of the Coulomb 
interaction disappears at T = in the case of the singlet 
ground state as discussed by Glazman and Matveevf^^ 
The Josephson current vanishes for the magnetic doublet 
ground state in the limit of |Afl|, |Al| — > cx). A finite 
current flowing in the opposite direction in the case of 
the finite SC gaps,'^ is caused by the continuous excited 
states above the gaps. 



At T = 0, the current can be expressed in the form 



eEAQ sin < 



2a/1 - Tsin^ 



Ea < 



U 



(Bll) 
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